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Abstract. An explanation is presented for the rather strong total surface differential rotation of the observed very young solar- 
type stars like AB Dor and PZ Tel. Due to its rapid rotation a nonuniform energy flux leaves the stellar core so that the outer 
convection zone is nonuniformly heated from below. Due to this 'gravity darkening' of the equator a meridional flow is created 
flowing equatorwards at the surface and thus accelerating the equatorial rotation. The effect linearly grows with the normalized 
pole-equator difference, e, of the heat-flux at the bottom of the convection zone. A rotation rate of about 9 h leads to e — 0.1 
for a solar-type star. In this case the resulting equator-pole differences of the angular velocity at the stellar surface, varies 
from unobservable 0.005 day -1 to the (desired) value of 0.03 day -1 when the dimensionless diffusivity factors c„ and c x vary 
between 1 and 0.1 (standard value c„ ~ c x ~ 0.3, see Table [j}) In all cases the related temperature differences between pole 
and equator at the surface are unobservably small. 

The (clockwise) meridional circulation which we obtain flows opposite to the (counterclockwise) circulation appearing as a 
byproduct in the A-theory of the nonuniform rotation in outer convection zones. The consequences of this situation for those 
dynamo theories of stellar activity are discussed which work with the meridional circulation as the dominant magnetic-advection 
effect in latitude to produce the solar-like form of the butterfly diagram. 

Key words. Hydrodynamics, Stars: rotation, Stars: pre-main sequence, Stellar activity 



1. Introduction 

In a series of papers the internal rotation of outer convec- 
tion zones of cool stars have been derived using a consis- 
tent mixing-length model of Reynolds stress and convective 
heat transport in a rotating convection zone that takes into ac- 
count the effect of the Coriolis force on the convective mo- 
tions (Kiiker et al. 1993, Kitchatinov & Rudiger 1995 (KR95), 
Kitchatinov & Rudiger 1999 (KR99) ). In the KR95 model, 
the deviation of the heat flux from spherical symmetry causes 
a small horizontal temperature gradient, which partly neutral- 
izes the rotational shear as the force that drives the meridional 
flow. As a result, the model reproduces the original Kiiker et al. 
(1993) result with only one free parameter, the mixing-length 
parameter omlt- 

KR99 found the total shear, 

Sfi = n cq — fjpoic, (i) 

to be only rather weakly dependent on the rotation period, O. 
Collier Cameron et al. (2001) have combined measurements of 
differential rotation for the pre-main-sequence stars RX J 1508- 
4423, AB Dor, and PZ Tel with the sample of main-sequence 
stars from Donahue etal. (1996) and argue that the discrepan- 
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cies between the KR99 result and the relation between shear 
and rotation rate derived by Donahue etal. (1996) may be due 
to the mix of G and K stars in the latter paper. They suggest 
that 5Q may indeed be constant with rotation rate but vary with 
spectral type and derive values for G and K dwarfs that dif- 
fer by a factor of three, the K dwarfs rotating more rigidly. His 
discussion of the sample of photometric data for several classes 
of stars has led Hall (1991) to the surprisingly flat distribution 
profile Sfl oc £1 015 . Further observations must decide how flat 
the 6Q — O relation for solar-type stars really is. 

In Fig. [I] the equator-pole differences ([!]) for the known 
rapid rotators are plotted in comparison to the results of the 
simulations by KR99. Obviously, the rapid rotators exhibit 
more differential rotation than the computations predict. One 
could believe - and we do - that another reason for differential 
rotation exists, which becomes important only for very rapid 
rotation. 

In the present paper the idea is adopted that for very fast 
rotation a nonuniform energy flux leaves the stellar radiative 
core ('gravity darkening'), thus heating from below the stellar 
convection zone nonuniformly (Stfpieft etal. 1997). This way, 
a meridional circulation should arise in the convection zone 
which even by itself produces a nonuniform surface rotation 
(Kippenhahn 1963). We shall see that an accelerated equator 
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Fig. 1. The total equator-pole difference of the surface angular 
velocity (1) versus the period r rot of the basic rotation after the 
model in KR99 and as observed 



is the natural consequence with an amplitude explaining in- 
deed the surprisingly steep surface rotation laws which are per- 
formed by the fast rotators. 



2. The Model 

In the diffusion approximation the radiative flux in a radiative 
stellar shell is proportionate to the gradient of the temperature 
which itself is a function of the total gravitational potential <j> — 

s 

ip — J n 2 s'ds'. Hence 



/(0)V0 = -f(<f>)g 



off 



err 



(2) 



Q 2 s. A con- 



with the effective gravity acceleration g 
vection zone above the radiative core is thus more intensely 
heated from below at the pole and it is less intensely heated 
from below at the equator. The intensity of this effect is given 
by the dimensional parameter 



Q 2 R 



(3) 



Its value for solar-type stars is e ~ 0.014/(r ro t/day) 2 with r ro t 
as the rotation period. For the very rapid rotation of 0. 1 days 
the e is of order unity and the star breaks off. In our computa- 
tions the rotation parameter e is considered as a free parameter, 
the models concern to a star with a rotation period of 9 h, i.e. 
e ~ 0.1. The same value would be true for a K giant with 19 
days rotation period. The KO giant KU Peg, however, has only 
e =0.04 (see Weber & Strassmeier, 2001, for differential rota- 
tion and meridional flown). 



1 Equatorial acceleration and poleward flow, quite similar as at the 
Sun 



2. 1. The transport equations 

We use the mean-field formulation of hydrodynamics, i.e. ap- 
ply an appropriate averaging procedure to split the velocity 
field into a mean and a fluctuating part. The mean velocity field, 
u, is then governed by the Reynolds equation, 



P 



^■ + («.V)« 



= -V • OQ) -Vp + pg, 



(4) 



where Qy = (u^u'j) is the correlation tensor of the velocity 
fluctuations, it'. 

As the density distribution is spherically symmetric, and the 
gas motion is dominated by the global rotation, we assume ax- 
isymmetry for the mean velocity and temperature. The velocity 
field can then be described as a superposition of a rotation and 
a meridional flow, u — r sin 9 fl cf) + u m , where cf> is the unit 
vector in the azimuthal direction. The azimuthal component of 
the Reynolds equation expresses the conservation of angular 
momentum: 



d (pr 2 sin 2 9 O) 
dt 



V-t = 0, 



(5) 



where t = pr sin (9 r sin 9 flu 111 + (u'^u 1 ) . As the mass den- 
sity varies with depth but not with time, mass conservation re- 
quires V • (pu) = 0. The meridional circulation can then be 
expressed by a stream function A. The stream function and the 
vorticity u> are related via the equation 



DA - 



IdpdA _ l_dp 
p dr dr pr dr 



A = 



-pu, 



(6) 



where D — A — l/(r 2 sin 2 9). 

The correlation tensor Q consists of a viscous and a non- 
viscous part. The eddy viscosity tensor, 



(7) 



has been calculated by Kitchatinov etal. (1994). It is here sim- 
plified to N = v^{8ik5ji + SiiSjk) ignoring the effect of the 
basic rotation to the tensorial structure. For the eddy viscosity 
jAp the mixing-length expression 



(8) 



u 2 ^ the turbulence intensity. The dimensionless parameter c„ is 
free but it should not exceed unity. Its standard value is 0.3 but 
under the influence of rotation also smaller values are possible. 
The non-viscous A-effect is the main source of differential ro- 
tation. In order to select the influence of the gravity darkening 
of the equator it is neglected, however, in the following calcu- 
lations. 

The convective heat transport is described by the transport 
equation 



ds 

pT— + pTu ■ Vs = -V • (F c 



F rad ) + q, 



(9) 



where q is the source function. In case of a perfectly adia- 
batic stratification the entropy is constant throughout the whole 
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convection zone. Standard mixing-length theory does not in- 
clude the rotational influence on the turbulent heat transport. 
Kitchatinov etal. (1994) have applied the same turbulence 
model as in their calculation of the Reynolds stress and found 



F = pTxtVs, 
where 



XT 



J/ 2 



(10) 



(11) 



is the mixing-length expression for the heat conductivity coef- 
ficient in a slowly rotating convection zone, i.e. fi* <C 1. Again 
c x is a free parameter. More details of the corresponding ther- 
modynamics are given by Kiiker & Stix (2001). 



2.2. Model of the convection zone and boundary 
conditions 

KR95 used a simplified model of the convection zone which is 
derived from the full model by solving the equations 

dT g(r) dg g T i 

from a starting point x e , where the reference values g e , p e and 
T e are taken from a standard solar model. Together with the 
opacity law, k — 0.34(1 + 6 x 10 24 pT~ 7//2 ) cgs, this gives 
a stratification quite close to that of the model the reference 
values were taken from. 

We additionally require that both the upper and lower 
boundaries are stress-free, 



}r<h — Qr6 = 0. 



(12) 



As boundary condition for the heat flux we require that the total 
flux through the outer boundary is equal to the total luminosity 



F r tot (r) = 



L 



47rr 2 



At the inner boundary the condition is now 



Airrf 



1 



e + 3 



(3 cos 2 9 - 1) 



(13) 



(14) 



The parameter e here gives the normalized difference between 
the heat flows at the poles and the equator, 



F r tot (pole) -F* ot (eq) 



(15) 



F* ot (eq) 

We have also tried the outer boundary condition of KR95 



F r tot (r e ) 



L (1 + 4 ™-) 



47rr 



CnTr 



pJ-cS 



(16) 



and found only small differences. However, ( jiff ) does not en- 
sure that the luminosity is constant with radius and the solu- 
tions show a slow drift in the total entropy while with (|l3]) a 
stationary state is reached after a few diffusion times. 
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Fig. 2. LEFT: The difference between the rotation rates at the 
equator and the the poles in deg/day. MIDDLE: The maximum 
value of the meridional flow velocity at the surface in m/s. 
RIGHT: The maximum value of the meridional flow velocity 
at the bottom in m/s. c v = c v = 1 



3. Results 

The results for c„ = c x = 1 are given in Fig. ||. We observe a 
nearly linear relation of the total equator-pole difference of the 
angular velocity and the amplitudes of the latitudinal drift at the 
bottom and the top of the convection zone with the darkening 
parameter e. We find 



U top - 31 £ 

and 

5tt ~ 0.058 e 



[day" 



(17) 



(18) 



Both effects, of course, are negligible for the present-day Sun 
(e ~ 10~ 5 ) but for rapid rotators with (say) e ~ 0.5 the effects 
should be observable. With (ph for solar-type stars it is 



Sn = 7.9-10 



-4 

'rot 



(19) 



For r ro t = 12 h one finds 5fl ~ 0.003 day -1 , a value which 
cannot explain the situation described in Fig. [l|. For rotation pe- 
riods exceeding 0.5 day and for the used turbulence model the 
gravity darkening effect seems not to play an important role. As 
we shall see below, remarkably higher values for the equator- 
pole-difference are possible for modified turbulence parameters 
but only up to a factor of 6. 

Another finding concerns the meridional flow. It has to do 
with the idea to interpret the solar butterfly diagram as the result 
of a meridional circulation (equatorwards at the bottom of the 
convection zone) rather than as an intrinsic property of the aVi- 
dynamo (Wang etal. 1991, Choudhuri etal. 1995, Dikpati & 
Charbonneau 1999). To this end the eddy diffusivity in the con- 
vection zone must be much smaller than currently believed. A 
typical value of a suitable magnetic Prandtl number with which 
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Fig. 3. The meridional circulation at mid-latitudes (45°) at the 
bottom (LEFT) and the top (RIGHT) of the convection zone for 
both a solar model and a K5 main-sequence star as a function 
of rotation period (KR99). The circle represents the present- 
day Sun. The flow at the bottom is equatorwards (as it is the 
flow at the surface, there are 2 cells in radius) 

Table 1. Meridional flow at top and bottom of the convection 
zone, equator-pole differences of angular velocity and tempera- 
ture for e = 0.1 (rotation period of ~ 9 h). Poleward circulation 
is marked with a minus 



c„ 


0.1 


0.1 


0.3 


1 


1 


c x 


0.1 


0.3 


0.3 


1 


0.3 


Pr 


1 


0.3 


1 


1 


3 


"top [m/s] 


+10 


+4.1 


+5.3 


+3.1 


+7.3 


Ubot [m/s] 


-1.3 


-0.46 


-0.81 


-0.58 


-1.50 


50, [day- 1 ] 


0.03 


0.013 


0.012 


0.005 


0.010 


5T [K] 


2 


1 


1 


0.5 


1 



meridional flows of amplitude of 10 m/s become magnetohy- 
drodynamically important is about 50 (Kuker etal. 2001). The 
equatorward bottom drifts for stars with various rotation rates 
are given in Fig. |[ It varies from 5 m/s for the present-day Sun 
to 11 m/s for a very young solar- type star. The bottom drift for 
the gravity darkening effect scales as 

u bot ~ -5.8 • e (20) 

in m/s. It flows opposite to the meridional flow induced by the 
A-effect, without nonuniform heating from below. As a result 
of the nonuniform heating the influence of the meridional flow 
onto the migration of the stellar toroidal field belts is reduced. 
A new and interesting question is whether the total meridional 
flow can change its sign by the influence of nonuniform heat- 
ing. The amplitude (|20|), however, seems to be too small to 
overcompensate 10 m/s even for the fastest rotation. 



We have thus to vary the free eddy diffusion parameters c„ 
and c x — c„/Pr. The results for a variety of parameters c v and 
Pr are presented in Table RL In all cases they are written for e = 
0.1, i.e. for a rotation period of about 9 h. There are no strong 
dependencies of the results on the input parameters c„ and Pr. 
We have given in Table [I] also the heat conductivity parameter 
c x = c v /Pr. Indeed, the equator-pole difference of the rotation 
rate has an inverse trend with c x . For c x = 1 we find only 50% 
of SCI resulting for c x = 0.3. For c x — 0.1 the factor of the 
corresponding Sft is 6, very close to the observations. 

On the other hand, there is a clear trend of the meridional 
circulation at the bottom of the convection zone with the tur- 
bulent Prandtl number. When the Prandtl number grows by a 
factor of 10 the bottom drift increases by a factor of 3. One can 
also find the reduction of the meridional flow for decreasing 
viscosity first described by Kohler (1969). His models demon- 
strate a distinct maximum of 10 m/s for an eddy viscosity of 
about 5 ■ 10 13 cm 2 /s. Density stratification, however, is not 
included. Models with density stratification require a faster 
meridional flow (factor 4 for a density contrast of 200) than 
models without density stratification (Riidiger 1989, p. 113). 
For an eddy viscosity of 5 • 10 12 in a density-stratified con- 
vection zone an equatorward circulation of 8 m/s at the surface 
produces a equator-pole difference of the angular velocity of 
SCI ~ 0.06 day" 1 . The data given in the second column of 
Table |l| confirm this well-known result (see also Kippenhahn 
1963). 

In all cases the latitudinal differences of the temperature are 
unobservably small. The maximal value for e = 1 is only 20 K. 
This finding is in accordance to the finding that the convection 
zone is able for an effective screening of all large-scale tem- 
perature differences at the bottom of the convection zone (Stix 
1981). 
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